SOLUTION OF THE CAUCHY PROBLEM FOR THE NAVIER - STOKES AND EULER 

EQUATIONS 

A. TSIONSKIY, M. TSIONSKIY * 

Abstract. Some known results regarding the Euler and Navier-Stokes equations were obtained by different authors. 
Existence and smoothness of the Navier-Stokes solutions in two dimensions have been known for a long time. Leray [2] showed 
that the Navier-Stokes equations in three space dimensions have a weak solution. Scheffer [5], [6] and Shnirelman [7] obtained 
weak solution of the Euler equations with compact support in spacetime. Caffarelli-Kohn-Nirenberg [8] improved Scheffer's 
results , and F.-H. Lin [9| simplified the proof of the results of J. Leray. Many problems and conjectures about the behavior of 
solutions of the Euler and Navier-Stokes equations are described in the book of Bertozzi and Majda [l] or Constantin [2]. 

Solutions of the Navier-Stokes and Euler equations with initial conditions (Cauchy problem) for two and three dimensions 
are obtained in the convergence series form by the iterative method using the Fourier and Laplace transforms in this paper. 
For several combinations of problem parameters numerical results were obtained and presented as graphs. 

1. The mathematical setup. 

The Navier-Stokes equations describe the motion of a fluid in {N = 2 or 3). We look for a viscous 
incompressible fluid filling all of i?^ here. The Navier-Stokes equations are then given by 

(1.1) ^ + Vwn|^ = '^Au, - ^ + fkix,t) G i?^, t>0, l<k<N) 

dt ^ dXn OXk 

n—1 



N 



dUn 



(1.2) divu = = {xeR'\t>0) 



n=l 

with initial conditions 

(1.3) u{x,0) = v°{x) {x e R^) 

Here u{x,t) = {uk{x,t)) e R^ , {1 < k < N) — is an unknown velocity vector (TV = 2 or 3), p (x, t) — 
is an unknown pressure, if'{x) is a given, C°° divergence-free vector field , fk{x,t) are components of a 
given, externally applied force f{x,t), is a positive coeflicient of the viscosity {ii v = then (1.1) - (1.3) 
are the Euler equations), and A = X^^^Li is the Laplacian in the space variables. Equation (jl.ip is 
Newton's law for a fiuid element subject. Equation (|1.2I) says that the fluid is incompressible. For physically 
reasonable solutions, we accept 

(1.4) Ukix,t) , ^ ^ as I X I oo (1 < k < N. 1 < n < N) 

OXn 

Hence, we will restrict attention to initial conditions ilP and force / that satisfy 

(1.5) \d"Tf{x)\<CaK{l+\x\y'^ on for any a and K > 0. 
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and 

(1.6) I 5^af /(x, t) I < C„0k(1+ I X I +t)-^ on x [0, oo) for any a, ^ and K > 0. 
We add (— J2n=i """ff^ ) *° both sides of the equations (|l.ll) . Then we have: 

(1.7) ^ = ^Az^fc - 1^ + A(x,t)- Vz.„|^ (xGi?^, i>0, 1<A:<7V) 

axfc oxn 

n—l 

N „ 

(1.8) divu = ^ pi = (x e > 0) 

n—l 

(1.9) u{x,0) = 77"(x) (x e i?^) 

(1.10) Ufc(a;,i)-J>0 , ^ ->Oas |x|->oo (l<k<N, l<n<N) 

OXn 

(1.11) I 9"it"(a;) I < CaK{l+ I x 1)"-^ on for any a and K > 0. 

(1.12) I d"d^f{x, t) I < Ca/3x(l+ I X I +t)"-f'^' on R^ x [0, oo) for any a, (3 and K > 0. 

We shall solve the system of equations (|1.7p - (|1.12p by the iterative method. To do so we write this 
system of equations in the following form: 

(1.13) ^ = v/^u,u - ^ + fok{x,t) (xei?^, i>0, l<k<N) 

(1.14) divu, = = (xei?^,t>0) 

n—l 

(1.15) itj(x,0) = if{x) (xeR^) 

(1.16) Ujkix,t)-^0 , -^Oas |x|-j>oo (l<k<N, l<n<N) 

(1.17) I 9"ti°(x) I < C„k(1+ I X \)-^ on R'^ for any a , K > and C^k > 0. 
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(1.18) I d^d^f{x,t) I < C„;3k(1+ I X I +t)-^ on x [0, oo) for any a,/3 , K > and C^fiK > 0. 
Here j is the number of the iterative process step (j — 1,2,3,...). 



N 



(1.19) M^,t) = Mx,t) - ^u,_i,„^|^ il<k<N) 

71—1 

or the vector form 



(1.20) fj{x,t) = /(x,t) - (uj_i • V)u,_i 

For the first step of the iterative process (j = 1) we have: 

(tio • V ) Mo 

and 

h{x,t) = /(x,t) 

2. Solution. Case N = 2. 

We use Fourier transform \K.2\ for equations (|1.13p — (|1.20p and get: 

f[ = -,ff/,..(,„,„.) [use«] 

(7r(7i,72) = FK(xi,X2)] 

Pj{jl,l2,t) = (xi,X2,t)] 
Fjk{ll,l2,t) = i^[/jfc(xi,X2,t)] 

fc,s = 1,2 



and then: 

(2_^) 9^j-i(7i,72,^) ^ -t.(72 + ^2)^^.^(^^^^^^^) z7iP,(7i,72,t) + F,i(7i,72,0 

(22) dU,2hu72,t) ^ _^(^2^^2)j;^.^(^^^^^^^) ^ Z72P,(7i,72,i) + ^,2(71, 72,0 
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(2.3) llUji{ji,-f2,t) + j2Uj2{jl,l2,t) = 

(2.4) C/ji(7i,72,0) = C/i"(7i,72) 

(2.5) t/,2(7i,72,0) = [/0(7i,72) 
Hence eliminate ^^(71,72,^) from equations (|2.ip . (|2.2I) and find: 

^ r 72 1 

7Jt[ ?7j2(7l,72,0 C/jl(7l,72,0J = 

at 7i 

(2.6) -K7f +7l)[ t/.2(7i,72,0 - ^;7,i(7i,72,i)] + [ ^^-2 (71 , 72 , - -i^,i(7i, 72, 

71 71 

(2.7) 7iC^ji(7i:72,0 + 72C^j2(7i,72,i) = 

(2.8) C/ji(7i,72,0) = C/?(7i,72) 

(2.9) C/j2(7i,72,0) = i7°(7i,72) 

We use Laplace transform (jA.4[) . (lA.sp for equations (|2.6p . (|2.7p and have: 

t^fft(7i,72,r/) = L[U,khul2,t)] k=l,2 

f/f2 (71,72, r/) - (71, 72,77)] - [(7,2(71,72,0) - ^C/,i(7i,72,0)] 

71 71 

(2.10) -i'(7i+7l)[C/f2(7i,72,^) - ^t/|?(7i,72,^)] + [ ^;l(7i, 72, ry) - ^ (71 , 72 , ^ 

(2.11) 7iC/fi(7i,72,^) + 72C/f2 (71, 72,^7) = 

(2.12) C/,i(7i,72,0) = C/i"(7i,72) 

(2.13) C/,2(7i,72,0) = C/2°(7i,72) 
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The solution of the system of equations (|2.10p — (|2.13p is: 



(2.14) C/fi(7i,72,r/) 



hl-Pjf (71,72,77) - 7i72-Fj|(7i,72,??) +7|t/i°(7i,72) - 71 72^/3° (71. 72)] 

ill + ii)[v + Hli + li)] 



(2.15) [71(71,72,7?) 



[7i-?vl(7i,72,77) -7172^71(71,72,77) +7iC72 (71, 72) -7172^/1(71,72)] 



(7? + 7|)[^ + K7?+7l)] 
Then we use the convolution formula (|A.6p and integral (jA.7p for (I2.14p , (|2.15p and obtain: 



(2.16) 



, ^^ [7|J;-1 (71, 72, t) -7172^^2 (71,72, t)] 

I7i,72,ij - e ^^,^^2) 

+ e -''(^?+^^')*C/?(7i,72) 



(2.17) 



[/,2(7i,72,0 = / e [7?^.-2(7i,72,r) -^7i72f.-i(7i,72,r)] _^ 

+ e -''(^'+^-')*[/2°(7i,72) 



^'j(7i,72,i) is obtained from (1^ or ((^ : 



(2.18) 



-P;(7i,72,i) = 7 



[71 -Fj 1 (71 , 72 , + 72 -^^2 (71 , 72 , 
(7?+7l) 



Use of the Fourier inversion formula (|A.2p and find: 



^ ^^(^?+-7;)ft-r) [7|-P;-i(7i,72,r) -7i72fj-2(7i,72,T)] ^_ 

(7?+7l) 



+ e -'^(^'+^^')*c/!'(7i,72) 



e d7id72 = 



00 /-CXD 



7l 



47r2 (7? + 7I) 7o 



00 /-OO 



e -K7i+72)(t-^) / / e *(^i'^i+*^''^)/ji(ii,i2,r)dxi(ii2(ir ■ 



— 00 ^ —00 
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-1 poo pOC P'' /" 



J — oc J — oc J — oc J — oc 



(2.19) 



Snifji) + SM2) + Biul) 



2 poo poo 
Uj2{xi,X2,t) = J J 



00 o —00 



-^(^?+^;)(t-r) bl ^^32 (71,72, t) - 7172^^1 (71, 72, t)] ^_ 

(7f+7l) ^ 



+ e -'^(t'+^=')*[/2°(7i,72) 



• e -'(^i7i+X272) ^^^^^2 + 



+ 



1 /-OO nOO 2 /"i /"OO /'OO 



• e -i(^iTi+X272) ^^^^^2 + 



-| poo poo poc poo 

/ Q -'^iil+'iDi / e '(*^Ti+^=T=)u^(ii,i2)rfiirfi2e -'(^i^^+^^^=)rf7id72 = 

J — 00 J — 00 t/— OO'/— 00 



(2.20) 



52l(/,l) + 522(/,-2) + B{ul) 



Pj{xi,X2,t) = — 



00 /■•■DO 



[71 ^ji (71 , 72 , + 72 ^"^-2 (71 , 72 , t)] 



00 J —00 



(7? + 7l) 



e d7id72 
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/oo poo poo poo 



^71^72 



poo /"OO poo poo 

^ I I 72 



OO •/ — oo 



(2.21) = + 52(/,2) 

So, the integrals (I^J^)) - (I^T^ exist by the restrictions prT7|) , pTTSI . 
Here ^nO, 5*120, 5*210, 522(), >5i(), 5*20 ^-''S the integral - operators. 

5l2() - 521 

We have for the vector Uj from the equations (|2.19p — (|2.20p : 

(2.22) Uj = S ■ + B{if) , 
where S is the matrix - operator: 

5ii 5i2 

521 522 

We put fj from equation (|1.20p into equation (|2.22p and have: 

Uj=S-if ~ ( 7lj_i • V)7lj_i) + B{if) = 
= f - {uj-i ■ V ) Uj-i + S(u°) = 

(2.23) = ui — 5 • {uj-i ■ V) Uj-i 

Here ui is the solution of the system of equations (11.131) — (|1.20l) with condition: 

V"n|^ =0 k= 1,2 

71 — 1 

For j = 1 formula (I2.22p can be written as follows: 

(2.24) ui = ~S ■ h + B{ifi) , /i(x,t) = fix,t) 

If t ^ then ui^ iP (look at integral-operators 5, B{) - integrals ([2TT91) , (|2?20l) '). 
For j = 2 we define from equation (|1.20p : 



(2.25) 



f2{x,t) = fl{x,t) - ( Wi • V ) Ml 
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We denote: 

(2.26) f; = (ui ■ V) ui 
and then we have: 

(2.27) f2{x,t) - Mx,t) -f; 

Then we get U2 from (12:22)) . ((2:241) : 

(2.28) U2 ■ f2 + B{v^) = 5 • (/i - ,/t) + B{if) = 7ii - 
Here we have: 

(2.29) ul ■ 

If t ^ then u^^Q (look at mtegral-operator S - integrals (j^l^ , ([2201)). 
Continue for j = 3. We define from equation (|1.20p : 

(2.30) /l(x,t) = h{x,t) - {U2 ■V)U2 
Here we have: 

(2.31) {U2 • V) U2 = ((ui - ui) - ui) = f; + f; 

We denote in ((23T1) : 

it = - ("1 • V) U2 - (U2 • V) Ui + {U2 ■ V) U2 

and then we have: 

(2.32) h{x,t) = h{x,t) -f; -f; 

Then we get from (1222]), 

(2.33) z?3 = 5 • + S(u°) = S ■ (/i - /t - it) + 5(^^) = - 
Here we denote: 

(2.34) u; ^ s . f; 

If t -> then ul ^ (look at integral-operator S - integrals (|2T^ , (f2:20|) ). 
For j = 4. We define from equation (|1.20p : 
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(2.35) Mx,t) = h{x,t) - • V)U3 
Here we have: 

(2.36) {u^ ■ V) W3 - ((W2 - ul) ■V){U2- ul) = f; + f; + n 
We denote in 

fl = - (m2 • V) ul - [ul • V) U2 + ("3 • ^^3 

and then we have: 

(2.37) h{x,t) = h{x,t) -f* - it 
Then we get U4 from (E^H), ^M), (12211), UMl- 

(2.38) U4 ^ S ■ U + B{tf) - S ■ (/i - f; - fS - fl) + B{v^) = - u*^ ~ ul - u\ 
Here we denote: 

(2.39) ul ^ ~S ■ fl 

If t — ^ then u| — > (look at integral-operator S - integrals (|2.19p , (|2.20p ). 
For arbitrary number j [j > 2). We define from equation (|1.20p : 

(2.40) f,{x,t) = fi{x,t) ~ (uj_i • V)u,_i 
Here we have: 

j 

(2.41) (w,_i • v)^r,_i =J2fi 

1=2 

and it follows: 



(2.42) f, = h -jZfi 



1=2 



Then we get Uj from ^TF^, (^TM)) : 



(2.43) u, = ~S ■ f, + B{u°) - 5 . (/i - 5^^*) + i?(S") = ^^1 - 

/=2 ;=2 
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Here we denote: 



(2.44) u* = S ■ /; (2 < I < j) 

If t ^ then ^ (look at integral-operator S - integrals ((2?T9| , (IS^ ). 

We consider the equations (|2.24p - (|2.44p and see that the series (|2.43p converge for j ^ oo 

with the conditions for the first step (j = 1) of the iterative process: 



, dXn 
n—1 



and conditions 



(2.45) CaK < ^ , CapK < ^ 

Here CaK and CajSK are received from (|1.17p , (jl.lSp . 
Hence, we receive from equation (I2.23P when j oo: 

(2.46) Uoc = ui - S ■ (uoo • V) Hoc 

Equation (|2.46p describes the converging iterative process. 
Then we have from formula ()2.2ip : 

(2.47) Poo = ^l(/ool) + 52(/oo2) 

Here /oo = (/ooi,/oo2) is received from formula (|2.42l) . 

On the other hand we can transform the original system of differential equations (jl.7p — (|1.9p to the 
equivalent system of integral equations by the scheme of iterative process (I2.22p , (|2.23p for vector u: 

(2.48) w = ui - § • (u • V) u, 

where ui is from formula (I2.24p . We compare the equations (|2.46p and (I2.48P and see that the iterative 
process (|2.46p converge to the solution of the system (|2.48p and hence to the solution of the differential 
equations (fLT)) - (fO)) with conditions ((2^ . 

In other words there exist smooth functions Poo(x, t), Uooi(x, t) (i = 1, 2) on x [0,oo) that 
satisfy (HI]), ([L2]), ([Ol) and 

(2.49) Poo, Uooi e C°°(R^ X [0,oo)), 



(2.50) / |uoo(x,t)|2dx< C 

for all t > 0. 
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3. Solution. Case N — 3. 

We use Fourier transform (jA.Sp for equations (|1.13p — (|1.20l) and get: 



^[ -Q^^a \ = -IsUjkili, 72,13, t) [use(|1.16p] 



C^fc(7i,72,73) = F[ul{xi,X2,X3)] 

Pj{ll,'l2,l3,t) = F[pj {xi,X2,X3,t)] 
Fjk{jl,'y2,j3,t) = F[fjk(xi,X2,X3,t)] 

k,s = 1,2,3 



and then: 



(3.1) ^^•^^'''^^^J^''^^'^'* = -^^(71 +72 + 73)C^ji(7i, 72,73, + «7i-Pi(7i, 72, 73, i) + -?^ji(7i, 72, 73, 



(3-2) ^^•^^^^^^^^'^'^'^'^ = -i^(7? +72 +7f)C^j2(7i,72,73,0 + «72^'i(7i, 72, 73, i) + -^^2 (71, 72, 73, 
(3-3) ^^^^^^^'J'^''^^''^^ ^ -i^(7i +72 +7f)C^j3(7i,72,73,0 + «73-Pi(7i, 72, 73, i) + -^^^3(71,72,73,0 



(3.4) 7if^ji(7i,72,73,0 + 72C^i2(7i,72,73,0 + 73 f^i3(7i, 72, 73, = 

(3.5) t/,i(7i,72,73,0) = C/i(7i,72,73) 

(3.6) C/,-2(7i,72,73,0) = C/2(7i,72,73) 

(3.7) C/,-3(7i,72,73,0) = 1/3(71, 72, 73) 
Hence eliminate Pj(7i, 72, 73, from equations p.ip — p.3p and find: 
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^[ C/,-2 (71, 72,73, i) - C/ji(7i,72,73,0] = -^^(7? + 72 + 73)[ ^^2(71, 72, 73, - 

(3.8) - — J7,-i(7i,72,73,t)] + [ -Fj2 (71,72,73,*) - — J^ji(7i, 72, 73, 

71 71 

4[ C^j3(7i,72,73,i) - — C^ji(7i,72,73,0] = -^^(7? + 7I + 7|)[ t^j3(7i, 72, 73, - 
at 71 

(3.9) - — C/ji(7i,72,73,0] + [ -F'j^ (71, 72,73, i) - — ^"^1(71,72,73,0] 

71 71 

(3.10) 7if^ji(7i,72,73,0 + 72^^^2(71,72,73,0 + 73 C^i3(7i,72,73,0 = 

(3.11) t^ji(7i,72,73,0) = C/i(7i,72,73) 

(3.12) «7,2(7i,72,73,0) = [/2°(7i,72,73) 

(3.13) C^j3(7i,72,73,0) = ?73(7i,72,73) 



We use Laplace transform (jA.4p . (|A.5P for equations p.8p — p.lOp and have: 

(7ffc(7i, 72,73,??) - i[f/jfc(7i,72,73,i)] k= 1,2,3 

^[ C{,®(7i,72,73,'7) - —t^^ (71, 72, 73,??)] - [ C^j2(7i,72,73,0) - — J7ji(7i, 72, 73, 0)] = 
71 71 

-^^(7? + 72 +7|)[ C^rt(7i, 72,73, T?) - — C^fi (71, 72,73, ??)] + 

7i 

(3.14) + [ ^;l(7i,72,73,^) - -^;'!(7i,72,73,^)] 

71 

^[ ^jl(7i,72,73,?7) - ^C^®(7i,72,73,'7)] " [ t/,-3(7i, 72, 73, 0) - ^?7ji(7i, 72, 73, 0)] = 

-K7i +72 +73)[ C/f3(7i,72,73,'7) - -C/,1(7i, 72, 73, ^y)] + 

71 

(3.15) +[j;|(7i, 72,73,^) - -F,!(7i,72,73,^)] 
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(3.16) 7iC^ji (71,72,73,??) + 72 t(j|(7i,72,73,?7) + 73 C/,|(7i, 72, 73, ??) = 

(3.17) C/ji(7i,72,73,0) = J7{'(7i,72,73) 

(3.18) C/,-2(7i,72,73,0) = t/2(7i,72,73) 

(3.19) (7^3(71, 72, 73,0) = 173(71,72,73) 

In the usual way the solution of the system of equations (|3.14l) — p.l6p with formulas p.l7p — p.l9p can 
be rewritten in the following form: 



55,, , _ [(7I +7|)-F/f (71, 72,73, >?) -7i72Jjl(7i,72,73,>?) - 7173-^^1 (71, 72, 73, ??) 

f^,i(7i,72,73,r?) - (7? + 7l+7|)[^ + K7?+7l+7l)] 



,0 o^^ , £^1 (71, 72, 73) 

^ ' [r? + K7? + 7l+7|)] 



X _ [(73 + 7? )-Fj| (71 , 72 , 73 , ??) - 7273-Fjl (71 , 72 , 73 , ??) " 7271 J^j^ (71 , 72 , 73 , ??)] 

ty,2(7l, 72,73,^) - (^2+^2+^2)[^ + ^(^2+^2+^2)J 



(3.21) 



^2'(7l,72,73) 

iv + Hif + il + 7I)] 



^, . _ [(7? +7l)-Fj|(7i,72,73,>?) -737iJjf(7i,72,73,>?) -7372fj|(7i,72,73,??)] 

(^,3(71,72,73,^) - (^2+^2+^2)[^ + ^(^2+^2+^2)] 



(3.22) 



1^3(71,72,73) 



[77 + 1/(72 + 71+72)] 

Then wc use the convolution formula (|A.6p and integral (|A.7p for p.20p — ()3.22p and obtain: 

t^ji(7i,72,73,i) = 

e -'/(7?+7'+7')(*-^) [(7I +7l)-P'3i (71, 72,73, t) - 7i72-Fj2(7i,72,73,T) -7173^3-3(71, 72, 73, t)] 
' (7i'+7l+7|) 

(3.23) + e f/0(7i, 72, 73) 
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J7i2(7i,72,73,i) 



(3.24) 



(-^2^-y;^^2^(^_^-) [(7I + 7f )Fj2(7i, 72, 73, t) - 7273J'j3(7i, 72, 73, t) - 727i-Fji(7i, 72, 73, t)] 

(7? + 7I + 7|) 

+ e -''(^'+''^'+^^')*C/2°(7i,72,73) 



1^^3(71,72,73,0 = 

g -^(-^,^+-y,^+^.;)(t-^) [(7? +7l)-Fj3(7i,72,73,T) - 737iJ'ji(7i,72,73,T) - 7372-^^2(71, 72, 73, r)] 
/o (7?+7l+7|) 

(3.25) + e -"(^I'+^i+^D* f/0(7i,72,73) 

-Pi(7i,72,73,0 is obtained from (01]) or ([331)] : 



(3.26) P,(7i,72,73,0 - « 



bi-P'j 1 (71 , 72 , 73 , + 72-F'i2 (71 , 72 , 73 , i) + 73-^,3 (71 , 72 , 73 , i)] 



(7?+7l+7l) 



Use of the Fourier inversion formula (|A.3P and find: 



Uji(a;i,a;2,X3,i) 



(27r)3/2 



00 pOO pOO 



00 ^ — 00 7 —00 



e -'^(-fHiU7i)(t-r) [(72 +73)-P^ji(7i, 72,73, r)] 

(7? + 7I+7I) 



_^(-y2^^2^^2)(t_^) [7i72-Fj2(7i,72,73,t) +7i73-Fj-3(7i,72,73,t)] ^ 

(7?+7l+7l) ^ 



+ e -(71+72+73)* [/0(7i,72,73) 



^ p p p ^(^2+^2) 

87r3 



1 

8^ 



1 

8^ 



00 J —00 J —00 



00 />QO pQO 



(7? + 7I + 7,1) 



e -''(7r+72+73)(«-^) 



g -»(^i7i+=>^272+a;373) (^^^^^2^73 = 
g i(5l7l+S272+£373) 



00 /'OO /'OO 



00 J —oc J —CO 



•/ji (Si, £2,^3,''') dxidx2dx^dT 



7172 



00 ^ —00 ^ —00 



00 /'OO nOQ 



(7? + 7I + 7,1) 



e -''(7j+72+73)(t-^) 



g -z(xi7l+a=272+a;373) d'yid'y2d-/3 - 
Q i(Sl7l+S272+£373) 



00 poo poo 



CO — 00 —CO 



■/j2(il,X2,i3,T) dxidX2dX3dT 



CO ./ — CO J —00 



7173 



(7? + 7I + 7I) 



e -''(7?+72+73)(t-^) 



g -,(a,i7i+:.272+:r373) djidj2dj3 " 

Q ■1(^171+^272+^373) 



00 /"OO pOO 



00 •/ —00 -/ —00 
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■fj3{xi,X2, X3, t) dXidX2dX3dT 



-1 fOC pOO fOC pOO fOC fOO 

_|_ / / / 6 ~''('''i+T2+T3)* / / / e *(*iTi+S272+S373) . 

J —oo J —oo J —oo .J—ooJ—ooJ — oo 



Ui{xi,X2, X3) dxidx2dxs 



-i(a;i7i +0:272 +X373) 



d'yid'y2dj3 = 



(3.27) 



^ll(/,l) + Si2(/,2) + 5l3(/,3) + i^K) 



2 pOO f-OO pOO nZ 



-!/(7 



?+7g+7g)(t-r) [(^/l + "l)^72(-l- -2- "3, t)] 



(7i +7I+73) 



-/' 

Jo 



_^(^2_^^2_^^2)(j_^) [7273-Pj3(7i. 72, 73, t) + 727i-Fji(7i, 72, 73, r) 



(71+72+73) 



+ e -'^(^i+^=+^3)t[/0(^^,^2,73) 



1 

8^ 



00 fOC fOO 



00 J —00 — CXD 



7271 



(7? + ll + 73) 



e -i/(7^+7^+7^)(t-T) 



g -i(xi7l+X272+X373) ^^^^^2^73 = 
g i(xi7l+X272+X373) 



00 /'OO /-OO 



00 —00 —00 



■/ji (ii, ^2, is, t) dx\dx2dx3dT 



g -i(xi7i+X272+X373) d-ixd-i2d'^3 + 



+ 



' ' (73 + 7?) 



87r3 .]_^ (72 + 7I + 7I) 



e ->'(7j+72+73)(*-r) 



/OO /"OO /"OO 
/ / 
-00 J —oc J —00 



g i(xi7l+X272+X373) 



■fj2ixi,X2, X3, t) dXidX2dX3dT 



g -i(xi7l+X272+X373)^^^^^2d73 



-j^ f-OC f-OC f-OO 

87r3 (72 



7273 



/OO /■• 00 /* 00 
/ / 
-00 J —00 - 



^+7l+7|) 

■fj3{xi,X2, X3,t) dxidX2dX3dT 



CX) —00 —00 



g i(xi71+X272+X373) . 

g -^(X171+X272+X373) djid-f2d-f3 + 



+ 



87r3 



OC /'OC /'OC 



00 —00 —00 



e -^(7^+72+73)* 



OC /'OC /"OC 



00 —00 —00 



g i(xi7l+X272+X373) _ 



U2{xi,X2,X3) dXidX2 dX3 



-»(xi7l+X272+X373) 



(^71^72^73 



(3.28) 



S2l{fjl) + 522(/i2) + 523(/,-3) + B{u^^) 



Uj3{xi,X2,X3,t) 



(27r)3/2 



OC /■• 00 00 



00 •/ —00 •/ —00 



■ r'g -uh?+^Uyi)(t-T) [(7i + 7l)-P;3(7i, 72, 73, t)] 
7o (7i + 7I + 7I) 
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.^^^2+^l+^2^^t-r) [7371 (71 , 72 , 73 , t) + 7372^ j2 (71 , 72 , 73 , t)] 



(7?+7|+7|) 



dr + 



+ e -^("'i+''=+''3)tf/0(^,,^2,73) 



g -i(a;i7i+X272+X373) (Ijidj2dj3 



2 /-CX^ /-CXD /'CXD 

87r3 (72 



7371 



/ e -^('ri+i2+'ri)(t-r) / / / 

,t/0 J — 00 J —00 J 



' + 7I + 7I) 

•/ji ^2, 2:3, t) dxidx2dx^dT 



i(5i7i +5272+^373) 



1 

8^ 



oc oc f 00 



00 */ —00 J —00 



7372 

(71+72+73) 



1 (7?+7l) 



I I I /-„,2 x ,,2 „2 



e -H-i7l+X272+u^373)^^^^^2d73 - 

I'T /'CO /'OC /'OO 

/ g -i^(7?+7|+7|)(*-T) / / / Q ^(^171+^272+5373) . 
vO J — 00 J — 00 J — 00 

fj2{ii,i2,S;3,T) dxidx2dx3dT 6 -'(^^^1+^^^=+^^^^) ^71^72^73 + 

/•OO /'OO pOO 

Q -^(7?+72+7|)(t-T) Q i(xi7i+X272+X373) . 

— CSO J —00 J —00 

Q -i(xi7i+X272+X373) ^71^72^73 + 



(7f + 72 + 73) 

■fj3{xi,X2, X3, t) dXidX2dX3dT 



+ 



87r3 



00 roc roc 



00 •/ —00 •/ —00 



Q -i^(7r+72+7|)i 



00 /'OO /'OO 



00 —00 —00 



g i(xi7l+X272+X373) _ 



U3{xi,X2,xs) dxidx2 dxs 



Q -i(a.i7i+X272+X373) djid-f2d-f3 



(3.29) 



= 53l(/,l) + 532(/,2) + 533(/i3) + B{ul) 



^ POO POO POO 



[71 ^ji (71 , 72 , 73 , + 72-Fj2 (71 , 72 , 73 , t)] 



+ 



+ 



^ /"OO /"OO /"CXD 

87r3 (72 



73-Fi3(7i;72,73,0 
(7?+7l + 7|) 

71 



"'+72+73) 

■fjl{xi,X2,X3, t) dxidX2dX3 



(7i + 7I + 73) 
g -i(xi7i+x272+x373) d'yid'y2d'y3 = 

/OO poo poo 
j / i(xi7l+X272+X373) . 

-00 J —00 J — 00 

g -i(xi7i+X272+X373) (;^^rf^2(i73 + 



8^ 



00 pOO pOO 



00 —00 —00 



(7i'+7|+7|) 



00 /'OO /'OO 



00 —00 —00 



i(xi71+X272+X373) 



■fj2{xi,X2, X3, t) dXidX2dX3 



g -i(xi7i+X272+X373) djidj2d^3 + 
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I 

8^ 



OO J — OO J — OO 



73 



ill + 7l + 7|) 



OO /'C30 /'OO 



OO ^ ~CXD ./ — OO 



71 +5272+^373) 



(3.30) 



■f]?,{xi,X2,X3, t) dxidX2dX3 
Slif.l) + S2iff2) + 53(/,3) 



(^71^72^73 



So, the integrals (IX?/)) - (1530)) exist by the restrictions prT7|) , pTTS)) . 

Here 5ii(), ^isO, ^laO, ^2i(), ^22 0, 523(), ^aiO, ^32(), ^33(), ^i(), ^2(), ^3() are the integral - oper- 
ators. 

5l2() = 521 

S13O = 531 



523() = 532() 

We have for the vector Uj from the equations p.27p — (j3.29p : 



(3.31) 



where S is the matrix - operator: 



Sii 


S12 


Sl3 


S21 


S22 


S23 


S31 


S32 


S33 



We put fj from equation (|1.20p into equation p.3ip and have: 

= 5 • ( / - ( uj^, ■ V)?r,_i) + B{tl") = 

= S ■ / - B- {uj_i ■ V ) + B{xf) = 

(3.32) — ui ~ S ■ {uj^i ■ V) Uj-i 

Here ui is the solution of the system of equations (I1.13P — (|1.20p with condition: 

3 



^z,„^ = k=l,2,3 

^ OXn 

For j = 1 formula (I3.3ip can be written as follows: 
(3.33) ui = ~S ■ fi + Biv^) , hix,t) = f{x,t) 
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If t then ui (look at integral-operators S,B{) - integrals ([5^ - ([5^ ). 

For j = 2 we define from equation (|1.20p : 



(3.34) f2{x,t) = h{x,t) - (ui ■ \7)u^ 
We denote: 

(3.35) f; = (ui ■ V) ui 
and then we have: 

(3.36) f2{x,t) = Mx,t) -f; 
Then we get U2 from ((33T|) . ((333| : 

(3.37) = 5 • /I + B(^ZO) - § • (/i - it) + S(^iO) = 7ii - 
Here we have: 

(3.38) ul = S ■ 

If t then -u^ (look at integral-operator § - integrals (P??7l) - ^^). 
Continue for j = 3. We define from equation (|1.20l) : 

(3.39) h{x,t) = h{x,t) - {U2 ■V)U2 
Here we have: 

(3.40) {U2 ■ v)u2 = {{ui - ^2) • V ) (^ii - u;) = f; + f; 

We denote in ((3^ : 

it = - (Ui ■ V)U*2 - {ul • V) Ui + {ul ■ V) U*2 

and then we have: 

(3.41) h{x,t) = /i(x,t) -/t - it 
Then we get U3 from (P3T|) . ([3331), (gSSl): 

(3.42) us = ~S ■ h + B{v^) = 5 ■ (/i - it - it) + B{v^) ^ ui ~ 
Here we denote: 
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(3.43) ^3 = ^ • /3* 

If t — > then (look at integral-operator S - integrals (|3.27p — (|3.29p ). 

For j = 4. We define from equation (|1.20p : 



(3.44) Mx,t) = h{x,t) - {us ■ \/)u3 
Here we have: 

(3.45) (m3 • V) w3 = (("2 - ui) ■ V) - u*) ^ f; + n + n 

We denote in (P^ : 

fl = - ■ y)u; - (^3 ■ V) U2 + ("3 • V) ui 

and then we have: 

(3.46) hix.t) = h{x,t) -f* - it 
Then we get U4 from (f33T|) . ^M), ([3:381) . ([QS]) : 

(3.47) ^ S ■ (/i - it - it - /T) + B{n^) = u, ~ u*, ~ u*, - 
Here we denote: 

(3.48) HI = S ■ n 

If t then ul ^ (look at integral-operator § - integrals (|3.27p - (|3.29p l. 
For arbitrary number j {j > 2). We define from equation (|1.20p : 

(3.49) f,{x,t) = fi{x,t) - (mj_i • V)uj-i 
Here we have: 



J 

(3.50) (^2,_i • V)S,_i = 5]it 

and it follows: 



(3.51) i- - fl - £it 
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Then we get uj from (P3T|) . 



(3.52) u, = ~S ■ f, + Biv^) ^ ~S . {f, - J2fn + Bin") = - E< 

1=2 1=2 

Here we denote: 

(3.53) =~S . fl {2 < I < j) 

If t ^ then ^ (look at integral-operator S - integrals (p?27l) , (|3^ ). 

We consider the equations p. 331) - p.53p and see that the series p. 521) converge for j ^ oo 

with the conditions for the first step (j — 1) of the iterative process: 

j^uor^^^O k = l,2,3 

n— 1 

and conditions 

(3.54) CaK< ^ , C„;3K< i. 

Here CaK and Ca^i^ are received from ()1.17p , (jl.lSp . 
Hence, we receive from equation p.32p when j ^ oo: 

(3.55) Uoo = - -S* • (Uoo • V) Uoc 

Equation p.55p describes the converging iterative process. 
Then we have from formula p.30p : 

(3.56) Poo = ^l(/ool) + 5'2(/oo2) + S^ifooi) 

Here /oo = (/ooi, /oo2, /oos) is received from formula (I3.5ip . 

On the other hand we can transform the original system of differential equations (jl.7p — (|1.9p to the 
equivalent system of integral equations by the scheme of iterative process p.3ip , (I3.32p for vector u: 

(3.57) u — ui — S ■ (u ■ V) u, 

where ui is from formula p.33p . We compare the equations p.55p and p.57p and see that the iterative 
process p.55p converge to the solution of the system p.57p and hence to the solution of the differential 
equations (|1.7p — (|1.9p with conditions p.54p . 

In other words there exist smooth functions Poo(x,t), Uooi(x, t) (i — 1, 2, 3) on x [0, cx)) 
that satisfy ([TI]), ([12]), ([O]) and 

(3.58) Poo, Uooi e C°°(R3 X [0, oo)), 
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(3.59) / |uoo(x,t)|2dx<C 

for all t > 0. 

In the following chapters 4 and 5 we describe in further details examples of the solutions for the Navier- 
Stocks and Euler problems with various applied forces and different values of the viscosity coefficient v. 

4. Example of the solution of the Cauchy problem for the Euler equations by the described 
iterative method with a particular applied force (N = 2). 

We will consider an example of the solution of the Cauchy problem for the Euler equations {v = 0) for 
N = 2 and with initial conditions: 

(4.1) u{x,Q) = iPix) = {xeR^) 

Hence, and from formulas (|2.19l) . (I2.20p for arbitrary step j of the iterative process, it follows: 

u,i{x,,X2,t) ^ — / , , / / / e ^(*^^i+^^^^)/,i(ii,i2,r)d£idi2dr. 

J -oo J -oo \ll + I2) Jo 



00 — 00 



47r2 J_^J_^ ill + ll) Jo 



00 -/ —00 



(4.2) 



I / 7172 



Uj2{xuX2,t) = --^ / / / / e ^(*^Ti+*^T^)/jl(il,52,T)dxidi2dT. 



00 ^ — 00 



-00 (7? +7I) Jo J-c 

e -'("1^1+"^''=) d7id72 



00 pOG ^,2 pt pOO nOO 

/ e '(*iTi+*^T^)/,2(il,52,T)dxidi2dT- 

00 J —00 



1 /"^ ry2 



47r2 (72 + 71) _/g 
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(4.3) 

We convert the Cartesian coordinates to the polar coordinates by formulas: 

Xi = r ■ cosifi ; X2 — r ■ siivp ; 71 = p • cos?A ; 72 = p • sin-0 ; xi — f ■ cosip ; X2 = f ■ sinip; 

and obtain from formulas (|4.2p . (|4.3p : 



471-^ ./o ^0 Jo Jo Jo 



-irpcos{ip — ip) 



pdpdijj 



(4.4) 



1 poo p27T pt poo p27r 

Uj2{r,ip,t) ^ -—^ / smiPcosiP / / e '^p''"''^^-'''^ fji{f,>fi,T) fdrd^pdr ■ 
47r Jq Jo Jq Jq Jq 



■ e -irpco^ii^-f) p^p^^ 



1 pOO pZTT pt pOO p2TT 

+ ^/ / COSV / / / e '^P^°''^^-''^fj2{f:^.r)fdfd^d7 



(4.5) 

We have the applied force fj for arbitrary step j of the iterative process: 

(4.6) /,?(?, ^, r) = /,,(0e ™^-^/,.(r) , /^^(f, ^, r) ^ 

or 
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(4.7) 



where fjf(f,Lp,T) , fj^{f^(p,T) — radial and tangential components of the applied force, 
nj - separate circumferential mode, Uj — 0,1,2,3,... 

We take the radial and tangential components of the applied force (|4.6p . (|4.7p with condition (|1.18p 
For the radial component of the applied force we use De Moivre's formulas (jA.SP and have: 

(4.8) /,2(f,^,r) = /,,(f)e "^-^szn^/,.(r) = ^/,,(f)(e ^("^-i)^ - e '("^■+i)'^)/,.(r) 



We put the applied force components (j4.8p in formulas (|4.4I) , (14.51) , change the order of integration and 



find: 



87r2 



^oo /'27r poo /'27r 

IQ JQ JQ JO 



Jo 



Jo 



-ir pcos{ip — ip) 



pdpdtp 



f]r{T)dT 



(4.9) 



87r2 



oo p2-k 



2tt 



simPcOS^P / fjf{f) / e ^fP'^'"^iv-i')(^Q ^(n,-l)<p_^g ^("^■+l)^)d(^fdf 

^0 Jo Jo 



+ i 



COS^lp / /ji=(f) / e i^P'^o'>if'-i'){Q iin,-l)v_Q i('h+m^dlffdf 



Jo 
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fjr{T)dT 



(4.10) 

We denote: 



(4.11) 

and from formulas (|ITTU|) it follows: 



"jt(0 = / f]TiT)dT 




J,jrl{r,(p,t) = Ujrl{r,ip)Ujt{t) 



(4.12) 
where 



OO /'27r pOO p27T 

"'0 Jo Jo 



■ e 



-irpcos{ip — <fi) 



pdpdip 



OO p2tt 



'o Jo Jo Jo 



(4.13) 



1 

8^ 



OO />27r 



JO 



ipcosip / fjf{P) / e «*=p™^('P-'A)(e »K-i)'?_^e '^"^+^')^)dipMr ■ 
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cos 



-irpcos{ip — ip) 



pdpdtp 



(4.14) 



Let us denote internal integrals in (|4.13p . (I4.14p as /+(f , p, ip) : 



(4.15) 



We have two integrals here. Plus (+) is for the first part of each integral (|4.13l) . (I4.14p and minus (-) is 
for the second part. 

We substitute for tp: 9 — - ^ , dO = dip and receive: 



g ifpcos9+i{nj~l)e^Q _|_ g i(nj+l)ip / g ifpcose+i{nj+l)0^g 

Then we use the Bessel function's integral representation (|A.9|) and have: 

(4.17) /+(r,p, V) = 2^i("^-i)e '("^-i)'^J„,-i(fp) + 2^i("^+i)e *("^+i)'^J„,+i(fp) 

Put /+(r, p, V') from ()4.17p in formulas (|4.13p . (|4.14p , change order of integration and obtain: 



Ujrl{r,ip) = 

(4.18) 



oo poo 



1 

8^' Jo Jo 



Mr) 



sin'^'ipl+{'r, p, Tp)^i sinijjcosijjl-if^ p, ip) 



e -'''p''°'^'f'-'^^di^fdf-pdp 



1 



oo poo 



-sin^cosipI+{r, p, ■0)+i cos^^I- (f, p, 'ip) 



e -'''P''°<'l'-fU^pfdfpdp 



u,^,{rM - I Mr) 

(4.19) 

Then we group parts in brackets of formulas (14. 18^ . (|4.19p , use De Moivre's formulas (jA.8[) and the 
Bessel function's properties. And we get: 



(4.20) ujri{r,ip) = - 



Mr) smi^e -'''-P'°'^^-'P^+'''^^ d^ Jn,{fp)d?dp 
Jo Jo Jo 



(4.21) Ujr2{r,(p) = 



2tt 



Jo 



fjf{f) / COSiP e ~^rpcos(^^^)+^n,^p J„^{Pp) dfdp 
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We substitute iox ip: 6 — ip - ip , dO = in the internal integrals of formulas (I4.20p . (|4.21l) , use 
De Moivre's formulas (jA.8[) and the Bessel function's integral representation (|A.9p and have from formulas 

6201), 6211): 



(4.22) u,ri{r,^) = "^'^ 



f]f{r) Jn,{rp) dfdp 



(4.23) Ujr2{r,ip) = 



e ^^J„,+,(rp) - e -^^Jn^.Arp) 



f]f{f) JnArp) dfdp 



Let us denote: 



(4.24) 



oo poo 



Rj,nj-l,r{r) = / f.jf{f) Jn^{fp)Jnj-l{rp) dfdp 

Jo Jo 



(4.25) 



oo /"OO 



10 JO 

Then we have from formulas (I4.22p . (|4.23p : 



/jfC?") Jn,{fp)Jnj+i{rp) dfdp 



(4.26) 



(4.27) 



Then if = it follows from 623), 623), 6211), 6^ that Ujri(r, f) = Ujr2(r, (^) 

Ml = 7/2 = 0. 

In the equations bellow we will consider Uj > 1. 

We change the order of integration in formulas (|4.24l) , (|4.25p and obtain: 



and hence 



(4.28) i?j- „,-i,r(r-) = / fjfir) Jn,{fp)Jn,-i{rp) dp df 

Jo Jo 



(4.29) 



Rj.7i,+i,r{r) = / fjf{f) / Jn,{fp)Jn,+i{rp) dp df 



Internal integrals in formulas (|4.28l) , (|4.29l) are established by the discontinuous integral of Weber and 
Schafheitlin (fOOt ^2^. Then we have from 623), 623): 



(4.30) 
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(4.31) Rj,n,+iAr) = ^ / r"^fjr{r)dr 







Now we integrate the solution (14. 4p . (|4.5p by the tangential component of the applied force (|4.7p for 
(rij > 1). Then we use De Moivre's formulas (jA.Sp and have: 

(4.32) 

Hence formulas (|4.32p are the components fji and fj2 from the tangential component of the applied 
force ()4.7p . while formulas (|4.8p are the components fji and fj2 from the radial component of the applied 
force 

Let us put in formulas and do the operations as we did in - (|i??7)) ( > 1). 

We consider that /jip(r) • /jriT) is restricted by condition (ll.lSp and get: 

/>oo />oo ^ 

(4.33) i?j,„^_i,^(r) = - / / {fj^{f) f)^ Jn^(fp)Ju,^i{rp) dfdp 

Jo Jo 

oo /-oo , 



(4-34) ^n.+i.-^CO = - / / {fj0{r) r)^ Jn,{fp)Jn,+i{rp) dPdp 

Jo Jo 

Here () ~ = Hence we have: 



(4.35) u,^iir,ip) = - ^[i?,,„^._i,^(r) 6 '("^"i)'^ + i?j,„,+i,^(r) e 'K+D^] 



(4.36) u,v^2(r,^) = i[i?,-„^._i,^(r) e '("^-i)'^ - i?j,„,+i,^(r) e ^n.+D^] 
We change the order of integration in formulas (|4.33p . (|4.34p and obtain: 

(4.37) i?j,„^._i,^(r) = - / {fjv{f)f)f\ Jn,{fp)Jn,-i{rp) dpdf 

Jo Jo 



(4.38) Rj,n,+i,v{r) = {fjAf)r)p Jn,{fp)Jn,+i{rp) dpdf 

Jo Jo 

Internal integrals in formulas (|4.37p , (|4.38p are established by the discontinuous integral of Weber and 
Schafheitlin (|Al0l) [12]. Then we have from (j437l) . g^Hl): 
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(4.39) 



R 



'i,nj-l.,ip 



(r) 



dr 



(4-40) %n.+i,^M ^ rifiAr) ri r~"^dr~ 

' ^ Jo 

We have obtain formulas (14.81) - (|4.40p for an arbitrary step j of the iterative process and the apphed 
forces (|46)) or (|4J)) . 

Now we investigate the first step (j = 1 , ni = n = 1, 2, 3, ...) of the iterative process with the particular 
radial component of the applied force fi{x,t) [ look at (|2.24p ]: 

Mr, ^, r) = /i,(f )e ^"^/i.(t) , /i^(f , ^, r) = 

(4.41) /i,(f)=F„f"+ie -^""^ , /i.(r)=e -^"^ 

Fn, o'n " Constants. < < oo , 1 < /i„ < oo , 1 < (T„ < oo. 

Let us put the applied force (|4.4ip in formulas (|4.30p . (|4.3ip . integrate and then we have: 

/•oo p fn+lp -fJ,„f pco p „n-l 

Ri,n-iAr) = / df = Fnr--' / e "^"^fdr = ^L_r(2,/x„r) 

(4.42) 



(4.43) 

r(a, a;), 7(a, x) are the incomplete gamma functions 
Hence, and from formulas (|4.26l) . (|4.27p it follows by j — 1: 

(4.44) ui.i(r,^) = |[i?i,„_i,,(r) 6 + i?i,„+i,.(r) e ^("+1)'^] 



(4.45) u,r2{r,^) = !^[i?i,„_i,,,(r) 6 ^(""^^^ - i?i,„+i,.(r) 6 

and from formula (14. lip : 



(4.46) 



Ult{t) = 



f hr{T)dT^ j e 

Jo Jo 



''dr = — 7(l,CT„t) 
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And we have the velocity ui [look at (|2.24l) ] from formulas (14.12^ : 

Uirl{r,ip,t) = Uirl{r,(p) Uit{t) 
(4.47) Uir2{r,(p,t) = Uir2{r, (p) Uit{t) 

We obtain the following equations by performing appropriate transformations: 



(4.48) uirir,^,t) = |[i?i^„_i^,(r) + i?i,„+i,,(r)] e "'^ ?xit(i) 



(4.49) ui^{r,if,t) = ^[i?i^„_i^,(r) - i?i,„+i,,(r)] e "'^ uu(t) 

uir{r, if, t), uiip(r, ip, t) are the radial and tangential components of the velocity ui. 
We have from formula (|4.46p : 

(4.50) hrnuu{t)=0 
Hence, and from formulas (|4.47p — (|4.49p : 

lim uiri{r,(p,t) ^ 0; lim Uir2ir, ^,t) = 0; 

(4.51) 

lim uir{r,(p,t) = 0; lim uiip{r,ip,t) = 0; 
In other words the velocity ui satisfies the initial conditions (|4.ip . 

We use the asymptotic properties of the incomplete gamma functions r(a, x), 7(0?, a;) and from formulas 
(|4.42p - (|4.45p , (|4.48p . (|4.49p we have the velocity ui satisfies conditions (|1.16p (for r ^ 00). 
Let us continue investigation for the second step (j = 2) of the iterative process. 
Find /2 ("T, <p, t) = {/li, 722} - the first correction of the particular radial applied force fi{x, t) (|4.41l) . 
We have for from formula p.26p : 

(4.52) /21 = U,rl + U,r2 ^ 



(4.53) /22 ^ ^1.1 ^ + U,r2 



where Miri, 1*1^2 are the components of ui and were taken from formulas (|4.47p . 
We have here: 
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duiri{r,p,t) 


duiri{r,p,t) 


dr 


1 duiri{r,p,t) 


dip 


dxi 


dr 


dxi 


dp 


dxi 


duiri{r,(p,t) _ 


duiri{r,p,t) 


dr 


1 duiri{r,p),t) 


dip 


8X2 


dr 


dx2 


dip 


dx2 


duir2{r,(p,t) _ 


duir2{r,p,t) 


dr 


1 duir2ir,p>,t) 


dp 


dxi 


dr 


dxi 


dip 


dxi 


duir2{r,(p,t) _ 


duir2ir,p,t) 


dr 


1 duir2{r,p,t) 


dp 


dx2 


dr 


dx2 


dp 


dx2 


dr dip 


sinp 


dr 


dp 


cosp 


— = cosp , — 

OXi OXi 


r 


dx2 


= sinp , - — ■ 

OX2 


r 



(4.54) 

Hence, we use formulas (|4.44l) — (|4.47p for uiri{r,p,t), uir2{r,p,t) and have from (14.54^ : 
+ * [(n - l)Ei,„_i,,(r) e ^("-1)^ + (n + l)i?i,„+i,,(r)e Hn+D'.] (_f!!^) U^^(i) 



+ z [(n-l)i?i,„_i,,(r)e ^("-i)'^ + (7i + l)i?i^„+i,,(r)e ^(^+^)^]^^ uu{t) 
+ * [(n - l)i?i,„_i,,(r) e ^("-1)^ - (n + l)i?i,„+i,,(r)e ^("+1)'^] (-£1!^) 1 uit(t) 



+ z [(n-l)i?i,„_i,,(r)e '("-1)'^ - (n + l)i?i,„+i,,(r)e ^ | ^^^(t) 



(4.55) 
where 
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(4.56) 

-x^ mo . _ 

dxi ' dx2 ' ' dx 



Let us put uir2, %^ from formulas (I071) . (1135]) in formulas (H351) 



for /21 , /22 • 
After completing appropriate operations we have: 

(4.57) f;iir,cp,t) = |J[r2,2„_i,.(r) e ^(^"^^^^ + T2,2„+i..(r) e ^^^^^(i) 

(4.58) f;2ir,V,t) - ^[T2,2„-i,.(r) 6 ^(^"-i)'^ - T2,2„+i,,.(r) 6 Wi)^]T2„(t) 
where 

'72,2n-l,r(?') = [Rl,n-l,r{r) + Rl,n+l,r{r)] Ri^n-l,ri^^ ^ ~ ^)Rl,n-l,r{r) |-^^^^_^^^(^^-) _ 

T2,2n+l,r{r) = [i?l,n- l,r (j') + Rl,n+l,r{r)] Ri^n+lA'^) ^ (» + l)Ri^n+l,rir) _ i?i^„+i_r (r)] 

(4.59) 

(4.60) T2nit) = uUt) 

We use formulas (|4.42l) . (I4.43P for i?i,„_i,r(r) , i?i,„+i^r(r) and (|4.56p for i?'i,„_i,r('") 1 ■R'i,n+i,r('') then 
do appropriate operations for T'2,2n-i,r('") , T'2,2n+i,r('") and get: 

7-2n-i 2fn — 1) 

T2.2«-i,r(r) - i^,? [ -^"'^r(2,/.„r) + 3 2n+4 7(2ri + 2,/i„r)r(2,^„r) - 

^ e -^"'■7(2n + 2,A.„r)] 



r2.2„+i,.(r) - [^^e -^"'-r(2,^„r) - 3 2^+4 7(2n + 2,/.„r) r(2,^„r) + 

+ ^r+2 e -'^"'■7(2n + 2,Ai„r)] 

(4.61) 

For radial and tangential /^^^ components of the first correction f2{r,(p,t) of the particular radial 
applied force we have: 
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(4.62) fUr,V,t) = ^[r2,2„-i,.(r) +r2,2„+i,,.(r)] 6 ^ T^nit) = ^r2,2„,r(r) 6 ' ^"^ T2„(t) 



■ 2 • 1 

in 771 

(4.63) /2%(r, V,t) = — [r2,2„-i,.(r) - r2,2n+i,.(r)] 6 ^ T2„(t) = _r2^2„,^(r) 6 ' r2„(0 
Here 

(4.64) T2,2„,.(r) = - i^,? , tn+4 7(2^ + 2, f^ar) r(2, /i„r) 

2j,2ri-l ^ 

T2,2„,v^(r) - - F2 [ -''"''r(2,^„r) - 3 2n+4 7(2n + 2,^„r)r(2,A^„r) + 

^T+2 e -^"'^7(2n + 2,Ai„r)] 



(4.65) 



We compare the particular radial applied force /i from ()4.4ip with the first correction /| from 
f (|4.62p — ()4.65p ) of this particular radial applied force, and we have: 

(4.66) I /I I « I /i 
with condition 

(4.67) Fn < - 

n 

After the first step of the iterative process (j = 1) we obtained the velocity ui [ see (|4.47l) ]. 
Now we will calculate ^2 - the first correction of the velocity ui. Solution of this problem has two stages. 
On the first stage we find the part of the first correction itj^, corresponding to the radial component of the 
first correction of applied force from (|4.62[) : 

(4.68) fUr, V. t) = ^T^anAr) 6 ^ ^"^ T^nit) , /2%(r, ^, t) ^ 

On the second stage we calculate the other part of the first correction 1^2^ , corresponding to the 
tangential component of the first correction of applied force /^^^ from ()4.63p : 

■ 9 

(4.69) f;,{r, ^,t)^0 , f;Ar, ^, t) = _T2,2„,^(r) e ' ^"'p T^nit) 



In other words 
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(4-70) U* = U*^+U*^, U*^ = {u2rl,W2r2}, = {"2>pl , ^^2^2} • 

First stage: we use formulas (|4.26p , (|4.27l) for components Ujri{r,(p), Ujr2{r,f) and formulas 
, (li3T]) forflw,,.,i .fr). i?j-„^+i,^(r) forj = 2and then formulas (gia) , (glMl) for A*^ (r, , T2,2„,rW- 
We do appropriate operations and have: 

(4.71) uU{r,^) = n[i?2,2„_i,,(r) 6 ^(^"-i)'^ + i?2,2„+i,r(r) e *(2n+i)^] 

(4.72) ul,^{r,^) = in[R2,2n-iAr) 6 ^^'""'^'^ - i?2,2„+i,r(r) 6 *(2"+i)'^] 
Here 

(4.73) -R2,2„-i,r(r) = ^2 / ^2;^ — 



-.2 /T 



77- 

(4.74) R2.2n+iAr) = ^^^5;^ r^"r2,2„..(r)df 







Second stage: we use formulas (|4.35|) , (|4.36p for components Ujtpi{r,(p), Uj^2{r,ip) and formulas 
Km . 610l)fori?j-„^._i,<^(r), R. r,,^^ Jr) for i = 2 and then formulas (I03l) . dUl]) for /2%(r, y^, i) , T2,2nAr)- 
We do appropriate operations and have: 

(4.75) u*,^,{r,^) = - ^[i?2,2n-i,^(r) 6 ^(^-D^ + i?2,2„+i,^(r) 6 ^(^n+D'P] 

(4.76) u^^^ir,^) = i[i?2,2„^i,^(r-) 6 '(^-D^ _ i?2,2„+i,^(r) 6 
Here: 

(4.77) i?2,2„-i,^(r) = ^ dr 

(4.78) i?2.2„+i,^(r) = -^2^:2^ r2"(T2,2„.^(f).f),df 
Then we have: 

u*2i{r,'p)=ul^i{r,(p)+u*^-^{r,ip) = 
= [nR2^2n-iAr) - ^i?2.2„-i,^(r)] 6 ^(^"-i)^ + [ni?2,2„+i.r(r) - ^i?2.2„+i,^(r)] 6 



(4.79) 



34 



A. TSIONSKIY, M. TSIONSKIY 



= ^ [nR2an^iAr) ~ |i?2.2„-i,^ W] 6 _ j [ni?2,2„+i..(r) - ^i?2,2n+i,^(r)] 6 

(4.80) 

From formula (|4.1ip for j = 2 and formula (I4.60p we have: 

(4.81) U2t{t) = f T2„(T)dT = ^[t-— 7(l,a„t) + -^7(l,2(7„t)] 
Hence, and from equation (|4.12p it follows: 

'^*2i{'r,'^,t) = ul^{r,ip) U2t{t) 

(4.82) ul^ (r, ^,t) ^ ul^ (r, ^) U2t (t) 
After completing appropriate operations we have: 

u;^{r,ip,t) = {[ni?2,2n-l,r(r) - ^R2.2n~lAr)] + ["^2,2n+l,r (?") " ^i?2,2«+l,v (r)] } 6 ^^"V ^^^(t) 

(4.83) 

ul^r.Lp^t) = i {[nR2,2n^l,r{r) - -R2.2n-l,^{r)] - [ni?2,2ri+l,r(7') " ■^R2,2n+l,^ir)]} 6 '^"'^ U2t(<) 

(4.84) 

Here U2r(r, 1^9, i), W2^(r, (/3,t) are the radial and tangential components of the first correction ut, of the 
velocity ui and 



2^2^2n-l 



i F Tt T 

[nR2,2n-lA'^) - 2-^2,2n-l,¥>(?')] = 



r(l,2^„r) r(2,2^„r) 



1=0 



22 

(2n + 2)i+i2'+V^ 



00 

i=0 



r(/ + 3,2^„r) 
(2n + 2),+i2'+VL 



(4.85) 



[ni?2, 



2nH-l,r 



(r) 



2,2n+l 



2,n2 r 



22j,2n+l 



n(2n + 1) 

94n-2 ,,4n+2 



7(4n, 2^„r) 
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+ 2n+2 7(2ri, ^^nr) r 6 ' " 

(4.86) 

From formulas (14.821) or (|4.83l) . (|4.84p with properties of W2t(i) - (|4.8ip it follows: 

, Urn ulj^{r,(p,t) ^0; Urn ul2{r,V,t) = 0; 



, , Urn ul^{r,ip,t) =0; Urn u* {r,ip,t) = 0; 

and we have the velocity U2 = ui — [ look at (|2.28p ] satisfying the initial conditions (j4.ip . We use the 
asymptotic properties of the incomplete gamma functions r(Q!,x), ^{a,x) and from formulas (j4.85p . (|4.86p 
we have: the first correction and therefore the velocity U2 satisfies conditions ()1.16p (for r oo). 

Let us compare the solution (|4.47p or ()4.48p . (|4.49p for ui of the first step of iterative process with the 
first correction (|4.82p or ()4.83p . (|4.84p for U2 , which is received on the second step of iterative process. We 
see that 



(4.89) I u* I << I iTi I 
with conditions 

n 

(4.90) t < a„ 

By continuing this iterative process we can obtain next parts itg , u| , ... of the converging series for u. 
For arbitrary step j of the iterative process we have by using formula (|2.43p : 

j 

(4.91) Uj = Ml - y^M* 

1=2 

and then: 

Urn Uj = u 

^ ' ^ j 00 

where u is the solution of the problem (11.11) — (|1.6I) for = 0. 

Below we provide numerical analysis of these results for the following values of problem's parameters: 
Circumferential modes n = 1, 2, 3, 4, 5. 
an = 10. 
< t < 10. 
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Results were obtained for functions ui - (|iTf|l or dUS]), gUl) ; u*^ - or with 

calculations of the incomplete gamma functions [llj . 



U2 = U1—U2 and is shown in FIG. 4.1 - 4.5. The vector field U2 at distances r = 1, 2, 3, 5, 7 is represented 
by the dotted curves in left diagrams. The comparison of | ui | (dashed plots) and | \ (solid plots) in plane 
Lp = [0, tt], at distances < r < 50 is represented in right diagrams. This comparison shows | U2 \ « \ ui \ 
and is corresponding to the conclusion (|4.89p . 
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5. Example of the solution of the Cauchy problem for the Navier - Stokes equations by 
the described iterative method with a particular applied force (N = 2). 

We will consider an example of the solution of the Cauchy problem for the Navier - Stokes equations for 
N = 2 and with initial conditions: 



(5.1) u{x,0) = u"ix) = (xei?2) 

Hence, and from formulas (|2.19|) . ()2.20p for arbitrary step j of the iterative process, it follows: 

-1 /"Oo />oo 2 rt f^oo /•oo 

Uj^{x^,X2,t) = — / / rr^^ / e / / e '^^'-^'+'^--<-^ f,i{x^,X2,T)dx^di2dT ■ 



00 J —00 



1 / / 7172 



4^2 (72 + 71) 



00 — CXD 



(5.2) 



u,2{xi,X2,t) = -— j j -^^^J^^ j^e -(^l+-<l)i^-) I I e ^(^^■'^+*=-^^)/,i(Si,i2,r)d£idi2dT. 



— 00 ^ — 00 



• e -^i^i^i+'^^T^) ^71^72 + 



+ / rT^^iHr / e -"^^i+^^'^^'-^M / e'^^'^^+^'^^^fj2iii,S:2,T)dxidx2dT- 



47r2 (72 + 71) 7o 



00 •/ —00 



(5.3) 

We convert the Cartesian coordinates to the polar coordinates by formulas: 

xi = r ■ cosifi ; X2 — r ■ shiip ; 71 = p • cos?A ; 72 = p • sin-0 ; xi — f ■ cosip ; X2 = f ■ sinip; 

and obtain from formulas (|5.2p . (|5.3p : 
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Uji{r,(p,t) = ^ J J sin'^tlj J e (*-^)y J e "'P''°'^^-''''>fji{f,<f,T)fdfdifdT' 



■t f-OC /■•27V f-t f-OO /■•27T 

j J sin^cosi; J e -"p'^^-"^ J J e **^'"=°'*(^-'^Vj2(r,<^,T)f(if#dr ■ 



• e -irpcosi^-v) pfipd^ 



(5.4) 



1 pOO pZTT pZ pOO pJi'K 

Uj2{r,(p,t) = J J sini^cosi; J e -''p^^*-'''> J J e '^^'"''^^-'^^ fji{f,<f,T)fdfd'fdT 

■ e -irpcos{^-p) pdpd^ + 

-t poo p2i\ pt pcyD p2r; 

+ J cos^i^ j e -"p'^'-"^ J J e '^p''°''^^-'>'^fj2{f,(p,T)fdfd(pdT 

■ e -irpcos(^-v) pdpd^ 

(5.5) 

We have the apphed force fj for arbitrary step j of the iterative process: 

(5.6) fjfif, ^, t) = fjfif, T)e '"^"^ , fj^if, <p,t)=0 
or 

(5.7) f,f{f,0,T)=O , /,^(f,^,r) = /,^(f,T)e "^'^ 

where fjf{f,(p,T) , fj^{f,(p,T) — radial and tangential components of the appHed force. 
rij - separate circumferential mode, rij = 0,1,2,3,... 
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We take the radial and tangential components of the applied force (|5.6p , (|5.7p with condition (|1.18p 
For the radial component of the applied force we use De Moivre's formulas (jA.Sp and have: 



(5.8) /j2(r,^,T) = fjf{f,T)e ™^-^szn^ = l/^.,(f,^)(e ^K-D^ - e K^.+Dv^) 

We put the applied force components (|5.8p in formulas (|5.4p . (|5.5I) and find: 



(5.9) 



1 /■°° r 



/•OO p2-7T 








Jo Jo 


.Jo Jo 



Jjry 



— K / / sinibcosib 



2Tr 



ijr2 



{r, ip, t) 



pOO p2TT 


pt p 


I I sinipcosip 


1 g -Up'it-T) 


Jo Jo 


Jo Jo 



27r 



8^ 



poo p27T 


pt p 


/ / cos^ip 


/ e / 


Jo Jo 


."'0 "'0 



2Tr 



g ifpcos(>p-V)(g iK-l)y_g iir^j+T-)v)d(frdfdT 



g -irpcos{4,-ip) p^p^^ 



(5.10) 
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Let us denote internal integrals in (I5.9p . (|5.10l) as /+(f , p, t/j) : 



(5.11) 



We have two integrals here. Plus (+) is for the first part of each integral (|5.9p . (|5.10p and minus (-) is 
for the second part. 

We substitute 9 for if : 9 — (p - ^ , d^ = d(^ and receive: 



g if pcosd+i{nj -1)6^0 g i(nj+l)V / g ifpcos9+i{nj + l)9^Q 

Then we use the Bessel function's integral representation (|A.9|) and have: 

(5.13) /+(r,p, V) = 2^i("^-i)e *("^-i)'''J„^_i(fp) + 2^^("^+l)e ^("^"+i)'^J„,+i(fp) 

Let us put I+{f,p,ip) from (|5.13p in formulas (|5.9p . (|5.10l) , change order of integration and obtain: 



Ujri{r,V,t) 



87r2 



"'0 "'0 
2Tr 







sin iljl^{f, p,%ji) — i sin4'cos4>I-{f , P, "0) 



(5.14) 



Ujr2{r,(p,t) 



oo /"i 



1 

c27r r 







Jo 

sinipcosipl^ (f, p,^) + i cos^tpl^ (f, p, -0) 



-irpcos(il) — ip) 



dipfdrdrpdp 



(5.15) 



Then we group parts in brackets of formulas (|5.14p . (|5.15p , use De Moivre's formulas (jA.Sp and the 
Bessel function's properties. And we get: 



(5.16) 



2tt 



Jo 



e / fjf{f,T) / sinV'e di; JnAfp) dMrdp 



(5.17) 
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We substitute 9 foi ^p: ~ ip - (p , d9 = dip in the internal integrals of formulas (|5.16p . (|5.17l) , use 
De Moivre's formulas (jA.8|) and the Bessel function's integral representation (IA.9I) and have from formulas 
(I536l),(l537l): 



(5.18) 



Jo 



e ■"^Jn.+.Crp) + e -^'^J„,_,(rp) 



fjf{r,T) Jn,{fp) dfdrdp 



(5.19) 



2 7o 7o 



e ''^J„,+,(rp) -e -'^J„,_,(rp) 



fjf{f,T) Jnj{fp) dfdrdp 



Let us denote: 



(5.20) i?j„._i,,(r,i) 



e -'^'''(*-")j„^.-i(rp) / /jK^t) J„^(rp) dfdrdp 

^0 



(5.21) i?,-„^+i,,(r,i) 



/O JO 

Then we have from formulas (|5.18p . (|5.19p 



e -''f"^'~^\jn,+i{^p) / /,f(r,r) J„^.(fp) dfdrdp 
Jo 



(5.22) 



(5.23) Ujr2{r,^,t) = ^[i?^,„^._i,,(r,<)e '("^"i)'^ - %„^.+i,,(r, t) 6 ^("^+1)'^] 



Then if = it follows from (IOT|) . that Ujvi fr. Lp. t) = Ujr2{r,f,t) = and 

hence ui = M2 = 0. 

In the equations bellow we will consider Uj > 1. 

Now we integrate the solution (|5.4I) , (15. 5p by the tangential component of the applied force (|5.7p ( for 
> !)• Then we use De Moivre's formulas (|A.8p and have: 



/,i(r~,^,r) = -fMr,r)e '^^^sin^ = _-/^.^(f, r)(e '("."D^ _e 



/,2(f,(^,r) = /,^(f,r)e ™>cos^ = i/,^(f,r)(e '(".-D^ + e 



(5.24) 
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Hence formulas (|5.24p are the components fji and fj2 from the tangential applied force (|5.7p . while 
formulas (j5.8p are the components fji and fj2 from the radial applied force (j5.6p . 

Let us put (|5.24p in formulas (j5.4p . (jS.Sp and do the operations as we did in (|5.9p — (|5.23p ( > 1). 
We consider that /j<p(r, r) is restricted by condition (|1.18l) and get: 

(5.25) Rj.n,-iAr,t) = - e -''p"^'-^\Jn,-i{rp) {f,^{f,T)-f)^Jr,^{rp)dfdTdp 



(5.26) %„,+i.^(r,t) = - / Q -"'^^'-^^JnM^P) {fj^{f:T)-f)^J^^{fp)dfdTdp , 

Jo Jo Jo 

Here () = Hence we have: 

(5.27) w,^i(r,¥.,i) = - ^[i?,,„,-i,^(r,i) 6 ^("^"i)'^ + i?,;„,+i^^(r, i) e '("^+1)^] 



(5.28) u,^2ir,^,t) = i[i?,-„,-i,^(r,t)e ^("^"i)'^ - i?,- „^.+i^^(r, i) e '('^^+1)^] 

We have obtain formulas (15.81) - (|5.28p for an arbitrary step j of the iterative process and the applied 
forces ^EM or (l57)) . 

Now we investigate the first step (j = 1 , ni = n = 1, 2, 3, ...) of the iterative process with the particular 
radial applied force /i(a:, t) [look at (|2.24p ]: 

/i,(f , ^, r) = /i,(f , T)e , /i^(f , r) = 

(5.29) /i^(f,T) = F„f"+ie -'^"'=Vir(r) 
FmPn - constants. , < < oo , 1 < < cxa. 

Let us put the particular radial applied force (|5.29p in formulas (j5.20p , (|5.21l) and for the internal 
integral we have by using formula (jA.lip . ^llj : 



(5.30) /(p,r) = /i,(f,r) J„(r» df = F„/i.(t) f^+^e J„(r» = 6 

Now we put /(p, r) from formula (I5.30p in formulas (|5.20l) , (I5.2ip , change the order of integration and 
have by using formula (|A.12p , JJj : 

Ri.u-iAr,t) = F„ / hr{T) / e ^""^ ' Jn-i{rp)dpdT = 

Jo Jo l^MrJ 



F„r 



,t/ir(T)$(n + 2,n + 2; 



(5-31) = u....... -^ ^'^^ 



2pl Jo [Afil,y{t - t) + ly 
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nt poo 

i?i,„+i,,(r,t) = F„ / /i,(t) / e 
Jo Jo 



(2M^) 



Jn+i{rp)dpdT 



(5.32) 



i+i /■* 



A,(r)$(n + 1, n + 2; ) 



4/^2 ^(t-T) + l 



2(n+l)7o [4AiXt - r) + 

Here <I>(a, c;a;) is a confluent hypergeometric function [10] . 

We substitute y for r: y = (^^y , dy = p-^-J^i^^^dr and receive 



(5.33) i?i,„_i,,(r,i) 



-1 



8m> 



/ir(y) • ■ ^{n + 2, n + 2; ~filr^y)dy 



p „n+l l-l 

(5.34) Ri,n+iAr,t) = , / /i,(y)-y"-i-$(n + l,n + 2;-M2^2y^^y 



8AiXn + l) 



Let us denote firiy) — and get: 



(5.35) Ri^n-i,r{r,t) = 



-1 /•! 



y" •$(7i + 2,n + 2;-^2^2y)dy 



(5.36) i?i,„+i,,(r,t) 



Fnr" 



-1 



i^(rj + 1) 



y"+^ ■<^>{n+l,n + 2:-p.lr^y)dy 



We use formula (|A.13p for integrand in the integral (|5.35p and formula (|A.14p for integrand in the integral 
(|5.36p 10 , integrate and then we have: 



Rl,n-l,r{r,t) 



n-l 



*(7i + 1, TT, + 2: -7t44xtt) 
<^>{n + 1, n + 2; -^^.^) ^...^t+.^n^l 



(5.37) 



-Rl,n+l,r(?', t) 



$(n+ l,n + 3;-/4r^) 



$(n+ l,n + 3; 



(^lilvt + 1) 



n+2 



(5.38) 



Hence, and from formulas (|5.22l) . (I5.23P it follows for j = 1: 



(5.39) 
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(5.40) u,r2ir,ip,t) = ^[i?i,„_i,,(r,Oe ^("^^^^ - i?i,„+i,,(r, 6 
We continue operations and get: 

(5.41) Uirir,V,t) - -[Ri,n-iAr,t) + i?i,„+i,.(r,i)] 6 "'^ 



(5.42) Mi^(r,¥',t) = —[Ri,n^i,rir,t) - i?i,„+i,,(r, t)] 6 "'^ 

uir{r, if, t), ui^{r, if, t) are the radial and tangential components of the velocity ui. 

We use the properties of the confluent hypergeometric function $(a, c; x) and have from formulas ()5.37p — 

lim uiri{r,(p,t) = 0; lira uir2('i','P-,t) = Q\ 

(5 43) * °' ^'^ °' 

lim uir{r,ip,t) = 0; lim ui^p^r, ip,t) ~ 0] 

Hence we have velocity ui satisfies the initial conditions (|5.ip . 

Then we use the asymptotic properties of the confluent hypergeometric function $(0, c; x) ilOj and from 
formulas (|5.37p — (|5.42l) we have velocity ui satisfies conditions (|1.16p ( for r 00). 
Let us continue investigation for the second step (j = 2) of the iterative process. 
Find /2 (?', ^, t) = {/21, /22} " the first correction of the particular radial applied force /i(x, t) (|5.29p . 
We have for from formula ([2:26| : 

(5.44) /21 = ""in y uir2 

0x1 0x2 

(5.45) /22 = + "lr2 

where wiri, uir2 are the components of iti and were taken from formulas (j5.39p . (|5.40p . 
We have here: 

duirijr, ip,t) _ duiri{r,ip,t) dr ^ duiri{r,(p,t) dip 
dxi dr dxi d(p dxi 

duirijr, ip,t) _ duiri{r,ip,t) dr ^ duiri{r,(p,t) dip 
dx2 dr dx2 dp dx2 

duir2{r, ip,t) _ duir2{r,(p,t) dr ^ duir2{r,(p,t) dip 
dxi dr dxi dp dxi 
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duir2{r,(p,t) duir2{r,(p,t) dr duir2{r,(p,t) dip 





dx2 


dr 


0X2 


dip 


dx2 


dr 


dip 


sinip 


dr 


dip 


cosp 


dxi 


= COSip , - — = 

axi 


r 


dx2 


smp , - — 
0x2 


r 



(5.46) 



Then we use formulas (I5.39p . (|5.40p for uiri{r,(p,t), uir2{i',ip,t) and have from (|5.46p : 



duiri{r,<p,t) _ njr, 

^("-'^'^ + (n + l)i?i,„+i^,(r,t)e 



dxi 

i [{n-l)Ri^n~iAr,t) e 
n 



duiri{r,(p,t) "frn' 
d^2 - 2lL^i-' 



i [{n - l)Ri,n 



2 ^ 

-iAr,t)e + (n + l)i?i.„+i. 



l9Mlr2( 



+ i 



[{n - l)i?i,„. 



i(n+i).^j ^ simp -^ 



^^f{K..-.,.,(.<)e 

_ (n + l)i?i„+i,,(r,t)e ^("+1)^1 £^ 



9X2 

i [(n-l)i?i,„_i,,(r,i) e 



(5.47) 
where 

-^l,n-l,r(^' 



ai?i,„_i..(r,t) ^«(^-i)^""'L^ 2 2^ ^(»+i>"+2;- (4";,+i) ; 

$(n+l,ri + 2;-^„r ) 



dr 

Fnr" 



A^llv{n + 2) 



8/i4i.(n + l) 



(4^>t+ l)"+i 



$(n + 2,n + 3;- 



{Afiliyt + l)"+2 



(5.48) 



-^l,n+l.r(^i ^) ^ 



dRi^n+l,r{r, t) 

dr 



F„r" 



SfiMn + 2) 



$(n + 1, n + 3; -7t44xtt) 
$(n + l,n + 3;-At„r-^) 
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n+2 



4z/(n + 2)(n + 3) 



$(n + 2,n + 4;-^2^2) 



$(n + 2,n + 4;-,j#- 



(4/^2 J/t + l)»+3 



(5.49) 



Let us put uiri, uir2, ^f, , -^Jf , froHi formulas ([0^1) . (lOO)) . ((CTf| in formulas 



(EUl), (S for /li , 
After computing appropriate operations we have: 



(5.50) 



(5.51) 
where 



f;2{r,V,t) = ^[T2,2„_i,.(r,i) e ^(2"-!'^ -r2,2„+i,.(r,t) 6 



T2,2n-i,rir,t) = [i?i,„_i,r(r, t) + i?i,„+i^r(r, t)] i?;^ „_i ,.(r, i) 



(n-1) 



-Rl,n-l,r(7-, [i?l,„_l,r(?', t) - i?l,„+l,r(f , i)] 



?2,2n+l,r(7',i) = [i?l,„_ l^r (j^, ^ ) + -Rl,n+l,r (j', *)] ^i,„+i,r (''i 



{n + 1) 



-Rl,n+l,r('',0[^l,«-l,r (?",*) " ^l,n+l,r(?', *)] 



(5.52) 

For radial and tangential components of the first correction f2{r,(p,t) of the particular radial 
applied force we have: 

? 2 

(5.53) /2V(r, i) = ^ [T2,2„-i,r(r, t) + T^^n+iAr. t)] e ^ = —T^^^^^r, t) e ^ ^"-^ 



(5.54) /2%(r,^,t) = ^[T2,2„-i,.(r,t) -T2,2„+i,.(r,t)] 6 ' = ^^2,2„,v^ 
where ( see (I5.52p ) 



T2,2„,v^(r,i)e '2"^^ 



T2,2n,r{r,t) = [i?i,„-. i,^ (r, t) +i?i^„+i,r(r,^)] [-Ri,„_i,r(r, t) + i?i^„+i_r(r, t)] 
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--[{n- l)i?i^„_i,r(?', t) + [n + l)i?i_„+i,r(r, t)] [i?i,„_i,r(r, t) - Ri^n+iA'^^ *)] 

T2,2n,^{r, t) = [i?l,„-l,r(?', t) + i?l,„+l,r(?', t)] [Ri^n-lA^, t) ' ■Rl,ri+l,r(^: *)] " 

-- [(n - l)i?i,„_i,r(r, t) - {n + l)Ri^n+i,rir, t)] [i?i,„_i,,.(r, i) - Ri,n+i,rir, t)] 



(5.55) 



We use formulas ([QT]) . ([OS)) for i?i,„_i,r(r, i) , i?i,„+i,r(?-, i) and ([5^ for i?i.„_i.r(r, i) , 

^i.n+i r(^i ^) then do the appropriate operations for T2^2n-i,r{r,t) , T2^2n+i,rir,t), using formula (IA.15p , 
and get: 



72,2«,r(?',0 = 



_p2 , ^in 



(5.56) 



$(n + l,n + 3;-/i2r2) 



1>(n+l,n + 3;-(3^g^) - 

(4/i2 iy< + l)«+2 



(5.57) 



2^2,2ri,ip(^j i) — 



$(n,n + 2;-/4r2) - 



(4/i>t+ l)"+i 



<l>(n + 2,n + 3;- . y 



F„ - n- r 



2n-l 



16Ai6i^2(„_^ X)2(n + 2) 



(4^>i + l)"+2 

_ _ $(„+!, „ + 2;-7T4i^^ 



$(n + l,n + 2;-^f,r^) - 



(4^.2 + 1) 



n+l 



$(n+l,n + 3;-M'7-') 



$(n+l,n + 3;-7T^'^' 



(4/i2 + 1) 



n+2 



By comparing particular radial applied force /i from (|5.29p with the first correction from ( (|5.53p 
(|5.57p ) of this particular radial applied force we have: 



(5.58) 

with condition 



/2 I « I /l 



(5.59) 



Fn < - 
n 
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After the first step of the iterative process (j = 1) we had the velocity ui - formulas (|5.4ip , (|5.42p . 
Now we will calculate U2 - the first correction of the velocity ul. Solution of this problem has two stages. 
On the first stage we find the part of the first correction U2r, corresponding to the first correction from 
formula (|5.53p of the applied force: 



(5.60) 



On the second stage we calculate the other part of the first correction u^^ , corresponding to the first 
correction from formula (|5.54p of the applied force: 



(5.61) 

In other words 



flir, ^,t)^0 , f^r, ^, t) = ^T2,2„,^(r, t) e * 



(5-62) U* = = {u2rl,W2r2}, = {"2^1' "2>p2}- 

First stage: After completing appropriate operations we have from formula (|5.56p : 



'2 ^2n-l 



72,2n,r(?',0 



1 



16Ai6t/2(n+ l)2(n + 2) 



+ 1, n + 2; -^i^r^) • $(n + 1, n + 3; -^J.lr'^) 



$(n + 1, n + 2; -M^r^) . $(n + 1, n + 3; -— 1^ 



(4/i2zyt + l)"+2 



(4Ai>t+l)' 



2 2 



2^2 



(4^>< + l)2»+3 



• $(n + l,n + 2; - 



lvt + l)> ^ ^ (4/.>i + l)^ 



2 ^2 



(5.63) 



We take formulas (|5.22p , (|5.23p for components Ujri{r,(p,t), Ujr2{r,ip,t) and formulas (|5.20p . (|5.2ip for 
i?j-„^._i,^(r,i), i?j-„^.+i,^(r,t) for j = 2 and then formulas (fH^HOll . (I5%1) for f^^{r,ip,t) , r2,2„,r(r,t) . We do 
appropriate operations and have: 



(5.64) 



4i(r,^,t) = n[i?2,2„-i,.(r,i)e ^('"-i)^ + i?2,2„+i,r(r, t) e ^(^"+1)^ 



(5.65) 



^2r2 



ir,ip,t) = zn[i?2,2„-i,.(r,0e '(^"-i)^ - i?2.2„+i,r(r, t) 6 '(^"+1)^ 



After changing the order of integration we receive: 



^2 pt poo poo 

(5.66) i?2,2„-i.r(r,t) = y T2,2nAr,r) J e -''p"^'-^'>J2n-i{rp) J2n{rp)dpdfdT 
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^2 ft foo 



(5.67) R2,2n+iAr,t) = ^ I I T2,2„,,.(f,T) / e --"P ^'-''>J2n+i{rp) J2n{rp)dpdrdT 



^0 



Second stage: After operations with T2,2n,ip{r, t) from formula ()5.57p and using formula (|A.15|) . we obtain: 

(rr ( ,\ \ ^(^2.2n.y(r,0 -r) / _l T i ,\ 

(T2,2n,ip[r,t) ■ r)^ = — = T2,2n,ip,vV',t) +T2,2n,<^,r(?',i) 



(5.68) 

We denote here 



2, 2n, 



-F„ ■ n ■ r 



,2n+l 



$ (n + 1 , n + 3: - 7x-t4xtt 



'i, ,2 2^ *(" + 2'^ + 3;-(4^f^; 
$(n + 2,n + 3;-/i2r2) - - 



^2 . ^2n+l 



(4^>i + l)"+2 
n + 2; - 



2 2 



(4^2 + l)"+i 



4>(n + 3,n + 4;-7:r44xTT 



{Aixlvt + 1) 



ri+3 



(5.69) 
and 



2 2 



_p2,^.^2n^i I $(n+l,n + 2;- 

<I>(n + l,n + 2;-/x„r^) v f-,. ^ ^ 



, 1 , o 2 2N ^(»+1." + 3;- (4m2Wi) . 



(4^>t + l)"+2 



(4^>< + l)»+i 

= 2n • T2^2n,r{r,t) 



(5.70) 



T2,2n,r{r,t) is taken from formulas (|5.56p and (|5.63p . 

Then we do several operations and have from formula (j5.69p : 



'-2,2n,ip,ip 



_p2 . ^ . ^2n+l 
8/i>2(„+l)(^ + 2)2 



$(n + 1, n + 3; -/i^'^^) • + 2, n + 3; -/i^^-^) 
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1 



(4^2^t + l)"+2 
1 

(4^>t + l)"+2 

1 



<I>(n+ l,n + 3;-/x„r^) • $(n + 2,n + 3; 



(4Ai>t + l)^ 



$(n+ l.n + 3; 



) •$(n + 2,n + 3;-/i^r2) + 



^(n +1 n + 3- ^-^i- 

(4/x2 + l)2n+4 + (4^2 + 1) 



(4m>< + 1) 

$(n + 2,n + 3 



2 2 



(4^>t+l)' 



^2 . ^2n+l 



Suti^^n + l){n + 3) 



$(n, n + 2; -/i^?-^) • *(?^ + 3, n + 4; 



2^2N 



"1^ 2 2 

- — — — r • + 2: -ulr^) ■ $(n + 3,n + 4: -- — 



1 



{A^ilvt + 1)^ 
1 

" (4^>i + l)2»+4 



— • $(n,ri + 2; 



(4/i2i.i+l)' 



$(n + 3,n + 4;-/i2 



2 2 2 2 

$fn.n + 2;- .. (f""^ . . ) ■ + 3, n + 4; 



(4M>i + l)^ 



(4M>i + l)' 



(5.71) 



We transform formulas (|5.27p . (|5.28l) for components Uj^i{r^Lp), Ujip2{r,(p) and formulas (j5.25p . (|5.26l) 
for Rj,nj~i.ip{f), Rj.rij+i.ipir) for j — 2 and then use formula (I5.6ip for /2<p(?', i)- We do several operations 
and have: 



(5.72) 



(5.73) u*2^2ir,^,t) = i[i?2.2„-i,^(r,t) 6 ^^^-D'. _ i?2,2„+i.^(r, t) 6 ^^^n+D'.] 



After changing the order of integration we obtain: 



2 /•* r°° 



(5.74) R2,2n-iAr,t) = - ^ / / (T2.2„,^(r,T) • f)^. / e J2„-i(r-p) J2„(rp) dpdfdr 



^0 



2 ft poo 



(5.75) R2,2n+i,v{r,t) = - — I / {T2,2nA^^^)-r)f I e ^'-^h2n+i{rp) J2n{fp) dpdfdr 



Then we take \T2^2n,ip{f,T) ■ f)^ from formula (j5.68p and with use of formula (I5.70p put it in formulas 
((5J4l) . ((5J5l) and have: 



in? /■* /■°° f"" 2 

R2,2n-l,y,{r,t) ^ ~ J J T2^2n,v,v{f,T) J Q '"^ ^*-~'''> J2n-l{r p) ■ J2n{f p) dpdfdr 
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in^ /■* f°° f°° 2, 

/ T2,2n,r(^, r) / G J2„-l (rp) • J2n (^p) (ipdflir = i?2,2n-l,<^,<^ - 2mi?2,2n-l,r 

^ Jo Jo Jo 
(5.76) 

f"" f"" 2 

R2,2n+l,v{r,t) = - — y y T2,2«,i^,i^(^,t) J G J2n+1 (j^p) " J2n{fp) dpdfdr - 



2 



T2.2n,r{r,T) 6 J2n+l{.rp) ■ J2nif p) dpdfdr = R2,2n+l,v,v{r,t) - 2niR2^2n+l,r{r, t) 

Jo Jo 



(5.77) 
where 



2 /■* /-oo 



(5.78) i?2,2n-i,^,^(^i) ^-^J^J^ T2,2„.^.^(r,r) 6 J2„-i(rp) • J2„(r» dpdfdr 

(5.79) R2,2n+i,vAr,t) ^ ~ — J J T2,2„,^,^(f,T)y G ""^ J2„+l (rp) • J2„ (fp) dpdfdr 

and i?2,2n-i,r(?', i) , ^2,2n+i,r(f, i) we take from formulas ()5.66p . (|5.67p . 
Then we use formulas ([QS]) . dSj!]), (fOSj) . (|5?72l) . ([STTS]) and have: 

"21(^7 i) = U2ri{r,(p,t) +U2^j^{r,(p,t) = 

= [nR2,2n-iAr,t) - ^R2,2n-iA'r,t)] e ^('""'''^ + [ni?2,2„+i,r(r,i) - ^i?2,2„+i,^(r,i)] e '(2"+!)'^ 
(5.80) 



■"22(''i <P) = W2r2(^: *) + ■"2^2(''i t) = 

= I [nR2.2n-iAr,t) - '-R2an-i..v(r,t)\ e _ , [ni?2,2„+i,r(r, t) - ^i?2.2„+i,^(r,i)] e 

(5.81) 

We obtain by performing appropriate transformations: 

U2r{r,(p,t) = { [ni?2.2n-l,r(f,i) " ^^2,2n- (?", t)] + [nR2an+l,r{r,t) - ^R2,2n+l,v{r,t)]} G '^""^ 



(5.82) 
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U*2^{r,(f,t) = i {[ni?2,2n-l,r - 2 -^2,2n-l,v ^)] " ["-R2,2n+l,r('', - 2 -^2,2n+l,<^ (f , i)] } G 

(5.83) 

Here M2r(f , "ys, ^), U2^{r^Lp,t) are the radial and tangential components of the first correction U2 of the 
velocity ui. R2,2n-i.r{r,t), R2,2n+i,r{T,t) are taken from formulas (|5.66p . (|5.67l) . 

R2,2n-l,v{r,t) = R2.2n^l,v,v{r,t) -2niR2^2n~l.r{r,t) 
R2,2n+l,vir,t) = R2,2n+l,ip,ip{.'r,t) ~2niR2,_2n+lA'^^^) 

(5.84) 

and R2,2n-i,.p,.pir,t), R2,2n+i,,p,,pir,t) are taken from formulas (|5.78p . (|5.79p . 
Then we do appropriate operations and have from formulas (j5.82p . (j5.83p : 

u*2^{r,ip,t) = -^[i?2,2«-i,v>,^(r,0 + i?2.2«+i,v,v(»-,0] e '^"^ 

(5.85) 

u*2Ar,ip,t) = ^[i?2,2„-i,^,^(r,i)-^2,2„+i,^,^(r,t)] e *2n<, 

(5.86) 

From formulas (|5.80p . (|5.8ip with properties of R2,2n~i,r{r, t), R2,2n+i,r{r, t), i?2,2n-i, ¥>(»■, *), -R2,2n+i,i^('', t) 
it follows: 



(5.87) 



(5.88) 



lim U2i{r,(p,t) = 0; lira U22{r,ip,t) = Q] 



lim U2r(r, Lp,t) = 0; lim W2^(r, 93, i) = 0; 
<-^0; t^O] 



In other words the velocity U2 = ui — [look at (12. 28^ ] satisfies the initial conditions (15.11) . We 
use the asymptotic properties of the confluent hypergeometric function $(a,c;a;) and have from formulas 
(j5.66p . ()5.67p . ()5.78p . ([5.79^ : the first correction and therefore the velocity U2 satisfies conditions (|1.16p 
( for r — 00). 

By comparing the solution ui from (|5.39p . (j5.40p or ()5.4ip . (|5.42p of the first step of iterative process with 
the first correction from (|5.80l) . (I5.8ip or (I5.85p . (|5.86p . which is received on the second step of iterative 
process, we see that 

(5.89) I u* I << I wi 
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with conditions 

n 

(5.90) 

By continuing this iterative process we can obtain next parts Ug ,^4..., of the converging series for u. 
For arbitrary step j of the iterative process we have by using formula (|2.43p : 

i 

(5.91) Uj — ui — W;* 

1=2 

and then: 



(5.92) 

where u is the solution of the problem (11.11) — (|1.6I) . 

Below we provide numerical analysis of these results for the following values of problem's parameters: 
Circumferential modes n = 1, 2, 3, 4, 5. 
/i„ = 1 (n = 1, 2, 3, 4, 5). 
Results were obtained for functions 

ui — (|5.39p . (|5.40p or (|5.41l) . (I5.42p with calculations of the confluent hypergeometric functions [lO]; 
u*2 - dEHOl), ([OT|) or (IES5)) . (jSTMl) by using numerical integration of the triple integrals ([E75)) . ([^Tg)) . 
Each of those integrals is computed as an iterated integral. 

Let us consider first the calculation of the inner integrals from (|5.78p . (|5.79l) : 

/•OO 

(5.93) I+{r,f,t,T)= e -•'p"^'-^^J2n+i{rp)-J2n{fp)dp 



For condition t < r the integrand is diminishing fast enough. It is easy to find upper limit of integration, 
so we can substitute integral (|5.93p for 

(5.94) I+{r,r,t,T)= [ 'e J2„+i(rp) • J2„(r» dp, 



where {0 < Ai = 200 < 00), and hence we are integrating over the finite interval. For additional check 
let us increase Ai in 1.5 times and change the number of steps of integration ni{ from 4001 to 6001 ). We 
see that the difference in result values of integral (|5.94p is within the range of required precision ei(10~"'^''^). 

For condition t = r the integral (|5.93|) is in fact an integral of Weber and Schafheitlin and it is possible 
to calculate it analytically (|A.10p 12 . 

Let us now consider the calculation of middle integrals 

/"OO 

(5.95) i+{r,t,T)^ / T2,2n,ifi,ifi{r,T)I+{r,f,t,T)df 

Jo 
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We use asymptotical properties of confluent hypergeometric functions $(a,c;x) '10[ and we have: 



(5.96) 



r ^ oo 



Hence, we substitute integral (|5.95p for 



(5.97) 



I+{r,t,T) 



{f,T)I+{r, r,t,T)df 



where (0 < A2 = 20 < 00) and integration is really over the finite interval. For additional check let us 
increase value A2 in 1.5 times and change the number of integration steps n2(from 201 to 301). We see that 
the difference in result values of integral (|5.97l) is within the range of required precision £2(10"^"'^). 

Confluent hypergeometric functions $(0, c;x) were computed with precision e(10^^^). 

The outer integrals in (|5.78p . (|5.79|) are the integrals over finite interval (0, t = 10). These integrals are 
computed with precision 63(10"^) and the number of steps of integration 77,3= 101. For additional check 
let us change the number of integration steps n3(from 101 to 201), and we see that the difference in result 
integral values is within the required precision 63(10^^). All integrals were computed by Simpson's method 
and ei{10-^^) < £2(10"") < esilO'^). 

U2 = ui — ul and is shown in FIG. 5.1.1 - 5.1.15. The vector field U2 at distances r = 1, 2, 3, 5, 7 is 
represented by the dotted curves in top diagrams. The comparison of | wi | (dashed plots) and | ul \ (solid 
plots) in plane (p = [0, tt], at distances < r < 50 is represented in bottom diagrams. This comparison 
shows I U2 I << I iTi I and is corresponding to the conclusion (|4.89p . 




FIG.5.1.1. n = 1, Fi = 1, ly ^ 1.5 
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FIG.5.1.5. n = 2, F2 = 0.5, u = 1 
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FIG.5.1.9. n = 3, ^3 = 0.33, v = 0.75 
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FIG.5.1.11. n = 4, ^4 = 0.25, v = 1 
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FIG.5.1.13. n = 5, i^s = 0.2, v = 1.5 
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FIG.5.1.15. n = 5, 



= 0.2, V = 0.75 
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Appendix A. 

The Fourier integral can be stated in the forms: 



N = 1 



-1 poo 1 poo 

(A.l) U{j) = F[u{x)] = j^-^ J^^ u{x) e ^-^^dx u{x) = f/(7)e -"'^dj 



N = 2 



1 pOO pOO 

«7(7i,72) = F[uixi,X2)] = — / u{xi,X2)e '(Ti^i+^^^^)rfa;i(ix2 

J — oo J — oo 
-I poo poo 

(A.2) u{xuX2) = — / [/(7i,72)e -*(^^^^+^^^^)d7irf72 



N = 3 



-j pOC pOO pOC 

U{jl,72,ls) = F[u{xuX2,Xs)] = j^;^ J J J u{x„X2,xs)e'^^'^'+'^'^'+'^'^'Uxidx2dxs 

-j^ pOO nOO pOO 

(A.3) u{xi,X2,X3) = ^2^^3/2 J J J f^(7i,72,73)e 



r-cjo poo poc 

-i(7ixi +720:2 +732:3) 



(i7id72(i73 



The Laplace integral is usually stated in the following form: 

/•oo -t /-c+ioo 

(A.4) u®{'n)=L[u{t)]^ u{t)e-''*dt u{t) = — u®{'n)e "^'dr] oco 

Jo 27ri Jc-ioo 

(A.5) L[w'(t)] =j?[/®(j?)-w(0) 

The CONVOLUTION theorem A.l. 
If integrals 

ufin) = / ux(t) e -"^'dt uf{ri) = / U2{t) e -"^dt 

Jo Jo 

absolutely converge by Rerj > ad, then U^{rj) — C/® (??) C/^f (??) is Laplace transform of 
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(A.6) U{t)= [ Ul{t -T)U2{T)dT 

Jo 

Useful Laplace integral: 

POO 1 

(A.7) L[e "'=*]=/ e = ^ [Rer, > %) 

Jo [n - Vk) 

De Moivre's formulas: 

(A.8) cosip = ^(e ^'^ + e -^^) , sirup = l(e - e -'^) 

Bessel function's integral representation: 

(A.9) Jn{z) = - — / e «cose+m0^^ ^ (j ^^-^ ^^^^^^ avouud the origin. 

27r Jc 

The discontinuous integral of Weber and Schafheitlin: 

(A.IO) / Jf,{at) ■ Jf,-i{bt) dt = \ 1/26 {b = a) 

■^0 I (6>a) 

(A.ll) / J^(at)e "'^ * rft= ^ i?e/x> -1, ^6 7^ >0. 

(A.12) Rej^>0,Re{ii + p)>0. 

(A.13) ^ [y° • $(a, c; -^j/)] = a ■ y''-' ■ $(a + 1, c; -/3t/) 

(A.14) ^ [y^-i • $(a, c; = (c - 1) • j/^-^ • $(a, c - 1; -^t/) 

Formula describing connection between the contiguous confluent hypergeometric functions: 

(A.15) c- $ - c- $(a-) - X • c • $(c+) = 
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